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A quantum inequality for the quantized electromagnetic field is developed for observers in 
static curved spacetimes. The quantum inequality derived is a generalized expression given 
by a mode function expansion of the four- vector potential, and the sampling function used 
to weight the energy integrals is left arbitrary up to the constraints that it be a positive, 
continuous function of unit area and that it decays at infinity. Examples of the quantum 
inequality are developed for Minkowski spacetime, Rindler spacetime and the Einstein closed 
universe. 
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I. INTRODUCTION 



Nearly four decades ago, it was shown by Epstein, Glaser and Jaffc jl] that a positive definite energy 
density was incompatible with the usual postulates of a quantized field theory. Worse yet, it appears that 
the energy density is not even bounded from below. Thus, all standard quantized field theories are capable of 
violating all the pointwise and averaged energy conditions in general relativity. However, this does not mean 
that the energy density can remain negative for an arbitrarily long period of time. Over the last decade, new 
forms of energy conditions involving various temporal and spatial averagings have been developed M pTJ . 
One such example is the quantum inequality, which is the weighted temporal average of the energy density 
along the worldline of an observer. Derived directly from quantum field theory, these inequalities limit the 
magnitude and temporal duration of existence of negative energy densities. The quantum inequalities say 
(3J), that if an observer tries to make a measurement of the energy density for some characteristic sampling time 

To, then the maximal negative energy that he might ever measure is bounded below by an inverse power of 
the characteristic sampling time. Given an observer's four-velocity u M and a sampling (weighting) function 
/ (t) of characteristic width tq, then the quantum inequality is given by 



/ a 
(Tfj,u(T)) Ren. U I ' l {t)u V (t) f{T)d,T > -S^To) + p va cuum- (1) 

-oo r 

Here a is a dimensionless constant of order unity and n is the dimension of the spacetime. For a massless 
field in Minkowski spacetime, the function S(tq) is equal to one and the vacuum energy density vanishes. For 
massive fields and/or curved spacetimes, S(t ) represents the modification of the quantum inequality away 
from its massless, flat space functional form. It has the generic behavior that it is approximately one for 
small To, and in most spacetimes it typically decays for longer characteristic sampling times. However there 
are some known exceptions, such as four-dimensional de Sitter and Rindler spacetimes, where the function 
S only grows only as fast as t§. 

The quantum inequalities were first derived by Ford |2] to constrain negative energy fluxes for the quan- 
tized, massless, minimally-coupled scalar field in Minkowski spacetime. These results were then expanded 
to the energy density of the massive scalar field in Minkowski space |3|,f7| and in static curved spacetimes 
P,^2p^]. In all of the these cases, a Lorentzian sampling function, 
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was used to simplify the calculations. However, Flanagan ]10| showed it was possible to derive optimum 
quantum inequalities for the massless scalar field in two dimensions for an arbitrary, smooth positive choice 
of the sampling function. This was followed by the work of Fewster and colleagues [^4 T5] T^] who have 
established the quantum inequality for the minimally coupled scalar field in static curved spacetimes of any 
dimension with an arbitrary, smooth positive sampling function. 

Although much of the previous work has been for the scalar field, work is now progressing for higher 
spin fields. Vollick has shown that an optimum quantum inequality can be derived for the Dirac field in 
two spacetime dimensions pfjfl for an arbitrary sampling function using the conformal properties of the field 
theory More recently, Fewster and Verch have established "quantum weak energy inequalities" for the Dirac 
and Majorana fields of nonzero mass in four-dimensional globally hyperbolic spacetimes Ell ]. Making use of 
microfocal analysis techniques, Fewster and collaborators P~yj2l| have vastly extended the applicability of 
the quantum inequalities to arbitrary globally hyperbolic spacetimes. 

The first quantum inequality for the electromagnetic field was derived by Ford and Roman J/J for a 
Lorentzian sampling function in flat spacetime. This was immediately generalized to curved static spacetimes 
by the author [Q , although both of these calculations relied on the specific choice of the Lorentzian sampling 
function. In addition, the proof in both cases was mathematically long and at some times quite complicated, 
particularly in some of the lemmas required. 

In this paper, we will show that it is possible to derive a generalized quantum inequality for the quantized 
electromagnetic field in static curved spacetimes with a length element of the form 



ds = — \goo(x)\dt + gij(x)dx l dx J . 



(3) 



The proof presented here is greatly simplified, in large part due to generalization of a more direct positivity 
lemma originally developed by Fewster and colleagues |14|,|l5| ■ In addition, the electromagnetic field quantum 
inequality is proven for an arbitrary choice of sampling function so long as it be a positive, continuous function 
of unit area that decays at infinity. The end result of our calculations is the quantum inequality written as 
a mode function expansion, 
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where Ej (A, k; x) and Bj (A, k; x) are the modes for the electric and magnetic components of the field-strength 
tensor, 



/ 1/2 H 
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is the Fourier transform of the square root of the sampling function and the summation over A and integration 
over d 3 k is over all possible polarizations and momentum eigenstates, respectively. As was the case for the 
scalar field, the electromagnetic field quantum inequality (0) tells us how much negative energy an observer 
may measure relative to the vacuum energy of the electromagnetic field. 

In Section O we will discuss the canonical quantization of the electromagnetic field in curved space and 
elucidate the particle state structure and the form of the stress-tensor. In particular, two different forms 
of quantization will be discussed: direct quantization in the classical Coulomb gauge and the more elegant 
Gupta-Bleuler form of quantization. In Section III we develop the positivity lemma for generic inner-products 
of vector fields, which is a generalization of work developed by Fewster and colleagues H,f5[ for the scalar 
field. In Section IV we lay out the remainder of the proof of the quantum inequality, finally arriving at 
the expression above. Lastly, in Section M we will look at the resulting quantum inequalities for Minkowski 
spacetime, Rindler spacetime and the Einstein closed universe. 

We will follow the the convention of Wald p2| where the signature of the metric is (—,+,+,+). Greek 
indices are summed over (0,1,2,3) while Latin indices denote the spatial components (1,2,3). However, the 



letter A has been singled out as the polarization state label, and depending on the context, can represent 
cither the two physical polarization states 1 and 2, or the full set of polarization states 0, 1, 2 and 3 in the 
Gupta-Bleuler formalism which includes the scalar and axial photon polarization states. Also, the complex 
conjugate of /, will be denoted by /. Units of fi = c = G = 1 will be used throughout. 

II. ELECTROMAGNETIC FIELD IN STATIC CURVED SPACETIMES 

We begin our discussion of the electromagnetic field by defining the classical Maxwell action for a source 
free field in curved space, 



^Maxwell _ 

4 Jv 



\J F aP F^^dd 4 x, (6) 



where F a p is the antisymmetric field-strength tensor related to the four- vector potential, A^, by 

F a p = X7 a Ap - y fj A a . (7) 

Here V represents covariant differentiation. 

Varying the Maxwell action with respect to the vector potential and setting the variation equal to zero 
leads to the source free inhomogeneous Maxwell equation for the electromagnetic field in curved spacetimc, 

V a F aP = 0. (8) 

Due to the Bianchi identities, the electromagnetic field also satisfies the subsidiary condition, 

V [a F 07] = O, (9) 

which is the homogeneous Maxwell equation. The combined set of equations represents classical electrody- 
namics written in covariant form. If we insert the four- vector potential into both expressions, it is found 
that the homogeneous Maxwell equation (||) is trivially satisfied. The inhomogeneous equation (gj) yields the 
second order wave equation 

V a V a A p - Vp [V a A a ) - Rp a A a = 0. (10) 

Here R a p is the Ricci tensor which arises due to the commutation relation for the covariant derivatives acting 
on a vector field. 

The stress-tensor for the classical electromagnetic field is found by varying the Maxwell action with respect 
to the spacetime metric. A straightforward calculation yields 

The field-strength tensor, the Maxwell equations and the stress tensor are invariant under the gauge 
freedom 

A ne W = A oU _ VqAj (12) 

where A = A(x) is an arbitrary scalar function. In classical electromagnetism, the correct choice of gauge 
can often simplify finding the solution to the field equations. In many cases, it is convenient to choose the 
Lorentz gauge condition 

which immediately removes the middle term in the wave equation (|I0|). This can always be achieved by 
choosing A to satisfy 

V"V Q A = V a A°J d . (14) 



It should be noted that there is still a restricted gauge freedom remaining in that we can still add to the 
vector potential any function that satisfies the homogeneous equation 

V a V a A Hom - = 0. (15) 

As we shall see below, this restricted gauge freedom will be used to impose the Coulomb gauge. It should 
be noted that we will drop the identifiers of new and old in all further calculations 

There is some difficulty is directly quantizing electrodynamics in the form so far described. If one does 
not specify a gauge, then any four-vector wave equation like (nG) will in general have four orthonormal 
solutions (polarization states), A a (X;x) where A = 0, 1, 2 or 3. In Minkowski spacetime the A = solution 
is typically the scalar photon polarization, A = 1 and 2 are the two transverse photon polarizations, and 
A = 3 is the axial photon polarization. In curved spacetime the "perfect" separation of the modes into these 
three "distinct" types is not always possible, but we will continue to use the flat space nomenclature. It 
is found that for one of the polarizations, say A a (0;x), there does not exist a conjugate momenta when 
the Hamiltonian is calculated. This is a long known problem in flat spacetime electrodynamics and there 
are several known approaches which have been developed to quantize the electromagnetic field that can 
be generalized to curved spacetime. The simplest is to work in a specific gauge [ p3|j24| . A more elegant 
possibility is to use the Gupta-Bleuler ]25W27J formalism of indefinite metrics on the Hilbert space of states. 
Both of these forms of quantization are discussed below. 

A. Direct quantization in the Coulomb gauge 

This is probably the simplest and most direct method of quantizing the electromagnetic field. The problem 
so far stems from the fact that the vector potential has four polarization states, while it is known that the 
photons of the free field theory only come in two different polarizations. Thus, before the theory is quantized 
we would like to remove the two superfluous polarizations at the classical level. To do this we require that 
solutions to the wave equation ( |l0| ) also satisfy the Lorentz gauge condition 

V a A a = 0. (16) 

This removes one degree of freedom between the components of the vector potential. The next condition 
that we would like to require is that the time component of the four-vector potential vanish in some frame. 
To accomplish this we let £ Q be a timelike vector field. Then we require that A a satisfy the additional 
condition 

CA a = 0. (17) 

It is this second condition that can be ensured by the homogeneous part of the gauge freedom. Also, note 
that is is not true that the Coulomb gauge is noncovariant as is sometimes stated. 

In flat spacetime there is no preferred choice of £ a , however for the static metric of the form (||) , a natural 
choice is to let £" to be the global timelike Killing vector field. This will be the same Killing vector that will 
be used to define the positive frequency mode functions. Since £ Q ex (1,0,0,0), the net effect is to set the 
Ao-componcnt of the stress tensor equal to zero. This solves two problems simultaneously. First it removes 
Aq from the action, thus there is no longer a problem of it not having a conjugate momenta. Secondly, it 
has reduced the physical degrees of freedom of the solution to the two physically realizable photon states. 

Canonical quantization is now straightforward. The metric (H) possesses a timelike killing vector, which 
allows us to write the the positive frequency mode function solutions of the wave equation ( |1C| ) as 

A Q (A,k;x,£) = £/ Q (A,k;x) e - *"*, (18) 

where k is the mode label for the propagation vector, A is the polarization state and oj = w(k). The four- 
vector functions, U^X, k; x), are the spatial portion of the solution of the wave equation and carry all the 
information about the curvature of the spacetime. In addition they satisfy 



V"[/ a (A,k;x) = = V Q £/ Q (A,k;x). (19) 



The mode functions for the vector potential are normalized such that 






(A(A,k),A(A',k')) = -i / d£ M A I/ (A,k)F^(A',k')-F^(A,k)^(A / ,k') - <5 AA 5 J (k - k'), (20) 



where cE^ = dan^ is a three- volume clement in the Cauchy surface £ with unit normal n M , thus each mode 
contributes \io to the vacuum expectation value of the stress tensor before renormalization. The general 
solution to the vector potential can then be expanded as 

2 . 

A ll (x,t) = J2 / d 3 k (a A (k) A„(X, k; x, t) + 4(k) A~{X, k; x, i)) . (21) 

Here at (k) and a A (k) are the creation and annihilation operators for the photon which obey the commutation 
relations 

[a A (k),a A ,(k')]=0=K(k),a|,(k')] (22) 

and 

[a A (k),4,(k')]=<5 AV( 5(k-k'). (23) 

The Fock representation of the number states can now be constructed from the vacuum state denoted by 
|0; 0) where the first slot is for particles of polarization type 1 and the second slot is for polarization type 2. 
The vacuum state has the property 

OA(k)|0;0)=0, V{A, k}. (24) 

One-particle states are obtained by acting on the vacuum with the creation operator, 

|l k ;0)=ot(k)|0;0) and |0; l k ) = 4(k)|0; 0). (25) 

Multi-particle states can likewise be created by repeated application of the creation operators, 

(aKki)) 1 " 1 • ■■( a t 1 (k 3 )) Jm ( a t(k 1 )) ln . . . (4(k,) 
( 1 m! . . . J m! 1 nl . . . in\) 



l 1 ™^, ..., 3 m^n^ 1 , ..., J n kj ) = ^ '- ^^ — " f - \ J ^ "-^- — 10;0), (26) 



where the ki , k2, . . . , kj are all distinct. The above state contains 1 m + 2 m + . . . + % + 1 n + 2 n + ...+ J n 
total particles where 1 m of them are of momentum ki and polarization 1, n arc of momentum ki and 
polarization 2, etc. Effectively, the general number states are a direct product of elements from two different 
Hilbert spaces, one for each of the polarization states. In order to reduce the index notation to a more 
manageable form, define the two vectors 

m= ( 1 mk 1 , •••, 3 rnk i ) and n = ^n^, . . . , J n^) , (27) 

then the states can be written more simply as |m; n). The most general state that can then be written as a 
linear superposition of all the possible number states, 

|V) = X)c(m,n)|m;n), (28) 

m.n 

where c(m, n) are complex coefficients and the sum is assumed to range over all the allowed vectors of m 
and n. For the state to be properly normalized, the c(m, n)'s must satisfy 

5>(m,n)| 2 = I. (29) 



qG.B. 



B. Gupta-Bleuler Formalism 

A more elegant form of quantization is to use the Gupta-Bleuler method of imposing an indefinite metric 
on the Hilbert space of allowable states. We begin by forming the Gupta action, 

nGupta nMaxwell , nG.B. f^fl^ 

where S Maxwel1 is the Maxwell action given by Eq. (||), and the gauge breaking action is given by 

l -J^ a A a f^-gd*x. (31) 

Variation of the new action with respect to A a yields the wave equation, 

V a F aP + V p (V a A a )=0, (32) 

which can be rewritten in terms of A a as 

V a V a A fj - Rf3 a A a = 0. (33) 

This would correspond to Maxwell's equations if the field also satisfied the Lorentz gauge condition. 

There are four possible solutions (polarizations) to the above wave equation. First, there are the two phys- 
ical polarizations which are labeled with A = 1 or 2. These two polarizations satisfy the wave equation ( |33| ) 
and the Lorentz condition, 

V Q A*(A,k;x,i) =0 for A =1,2. (34) 

Thus, these two polarizations correspond to the two standard solutions to Maxwell's equations. The re- 
maining two unphysical polarizations, labeled with A — or 3, also satisfy the wave equation (|33|), but not 
necessarily the Lorentz condition. For ultra-static spacetimes, where |<7oo| = L the most natural choice is to 
use the scalar photon polarization, 

A Q (O,k;x,t) = -(9 t ,O,O,O)0(k;x,t), (35) 

and the longitudinal photon polarization, 

A a (3, k; x, t) = — (0, d x i , d x 2 , d x a ) 0(k; x, t) , (36) 

U! 

where 0(k; x, t) is the solution to the massless, minimally couples scalar wave equation, 

V Q V Q (/.(k;x,t) = 0. (37) 

In the more general case of a static spacetime, it is useful to choose the two orthogonal modes which satisfy 
the condition, 

A a (0, k; x, t) + A a (3, k; x, t) = - V« <£(k; x, t). (38) 

In both cases, the resulting modes then satisfy 

V Q ^ Q (0,k;x,t) = -V a A Q (3,k;x,t) (39) 

and 

F Q/3 (0,k;x,i) = -F Q/3 (3,k;x,i), (40) 

for every momenta k. 

In addition, if we define the generalized conjugate momenta, 



IP" = -{F» v + g^V a A a ), (41) 

then the modes are required to be orthogonal and normalized by 

(A(\,k),A(\',k , )) = i /ds4^(A,k)n^(A',k')-n^(A,k)^(A',k')l = V xx 's 3 (k-k'), (42) 



diag(-l, 1,1,1). The general solution to Ay then has the Fourier mode-decomposition 

3 . 
A lt {x,t) = J2 / d 3 k (a A (k) Ay(X, k; x, t) + 4(k) A M (A, k; x, *)) . (43) 

If we wish to canonically quantize the field Ay,, we impose the equal-time commutation relations 

[Ay (x, t) , A v (x', t)\ = = [H> (x, t), n*" (x' , £)] (44) 

and 

[^(x, t), n*"(x', t)] = 4^-5 3 (x - x'). (45) 

Using the mode decomposition and the normalization condition, we find that the above equal-time commu- 
tation relations are equivalent to 

a A (k),4,(k')] =7 ?AA '£ 3 (k-k'), (46) 

with all other commutators vanishing. 

The state structure is similar in form to that found for the Coulomb gauge, except there are now a greater 
number of allowable states due to the two unphysical unphysical polarizations. We now define the vacuum 
state as |0; 0; 0; 0} where the first slot is for photons of the unphysical polarization A = 0, the second and 
third slots are for the two real photon polarizations, and the final slot is for the unphysical polarization with 
A = 3. The vacuum state vanishes if any of the four destruction operators act on it, and multi-particle states 
are again obtained by the repeated application of the creation operators. Unlike, the states for the Coulomb 
gauge quantization, the states of the Gupta-Blculcr formalism have indefinite norm, 

(l,m,n,p|l',m',n',p') = (-l) li+2 ' + - + ^n<S mm <(W<W, (47) 

where we have added two new vectors, 1 and p, for the unphysical photon polarization states. The most 
general state in the Gupta-Bleuler formulation can be written as superposition of all the particle number 
states as 

\(f>)= J2 c(l,m,n,p)|l,m,n,p). (48) 

l,m,n,p 

In order for the Gupta-Blculcr formalism to be equivalent to Maxwell's theory, we need to impose an 
additional condition on the Hilbert space of states; that the expectation value of the Lorentz condition be 
satisfied for all physically realizable states \<j>), 

(0|VM Q (x,i)|0)=O. (49) 

This condition can be accomplished simply by requiring that the states obey 

V a A+(x,t)\4>)=0, (50) 

where A^ is the positive frequency part of A a . The application of this condition to the state \<j>) above 
means that the c(l, m, n, p)'s with the same total number of A = and 3 photons of the same momenta are 
related to one and other by 



^VM^O, k)c (/ k , m, n, [p - l) k ) + VPk"V a A*(3, k)c ((I - l) k , m, n,p k ) = 0. 

Under this constraint the Hilbert space structure of the state \cf>) takes the form, 

|0) = hc(0 k ,m,n,0 k )|0 k ,m,n, k ) H 

c(l k ,m,n,0 k ) 



(51) 



V Q A Q (3,k) 

c(l k ,m,n, l k ) 

V2V Q A Q (0,k)V Q A Q (3,k) 



[VVL a (3,k)|l k ,m,n,0 k )-V Q A Q (0,k)|0 k ,m,n,l k )] + --- 

(V Q A a (3, k)) 2 |2 k , m, n, k ) - V2V a A a {0, k)V Q A Q (3, k)|l k , m, n, l k ) 

(52) 



-(V Q A Q (0,k)) z |0 k ,m,n,2 k 
With the definition of a new operator 

O f (k) = W a A a (3, k) aj(k) - W a A a (0, k) a|(k 



(53) 



it is possible to define a new set of states, |m, n, q), where each *g k . in q is the total number of A = and 3 
photons of momenta k^. The new states are formed by the repeated action of the operator O^(k) acting on 
the state with zero unphysical photons, 



I m > n >{ 1 9k 1 > 2 9k 2 ,---Akj) 
The inner product of the new states are 



(Ot(k 1 ))'(Qt(k 2 )) g ...(Qt(k i )) 
(ig! 2 ^!... * ff !)Va 



|0,m,n,0) 



(m, n, 0|m, n, 0) = 1 

for the states with no unphysical photons and 

<m, n, g k |m, n, q k ) = {\V a A a (0, k)| 2 - \V a A a {?>, k)| 2 } 9 = 



(54) 



(55) 



(56) 



for all other states. It is now possible to rewrite the superposition of particle number states ( ^8[ ) with the 
embodiment of the supplementary condition (|50| ) built in as 



10) = Yl fo ( m > n > c 0i m > n ' c i)- 

m,n,q 

The stress tensor found from the Gupta action is 

rpGwpta __ rpMaxwell ■ rpG.B. 
pa pa * pa 5 

-i Maxwell 



(57) 



(58) 



where T^ axwel1 is given by Eq. (pd[) and the contribution to the stress-tensor from the gauge breaking term 



t g.b. 



-A p (V CT V Q A Q ) - A a (V p V a A a ) + g pa 



\:\ i \7 a A a + ^(\7 a A a ) 2 



M' 



(59) 



Due to the physical photon polarizations modes satisfying the Lorentz condition ( |34| ) and the Hilbert space of 
states satisfying the the subsidiary condition (|o|) , it is relatively straightforward to show that the expectation 
value of the normal ordered gauge-breaking portion of the stress-tensor vanishes, 



:T,' 



G.B. 

pa 



= 0. 



(60) 



In addition, due to the relationships between the c(l, m, n, p) coefficients and Eq. (|4(j), it is simple to show for 
the normal-ordered Maxwell portion of the stress-tensor that the unphysical photon modes do not contribute, 
thus 



(0| : T£«»° M : |0) = (VI : T%°™ eU : \ij>), (61) 

where 

iv> = J2 c (°' m ' n ' °) i°' m > n ' °> = £ c ( m ' n ) i m ' n ) • ( 62 ) 

m,n m,n 

Thus, the only physically observable states are the two physical photon polarization states. In summary we 
have 

(0| : T^ ta : |0) = (0| : T%°*» M + Tf a B - : |0) = (0| : Tp A ^ e " : |^>. (63) 



III. POSITIVITY RESULT 

In this section we prove the following inequality: Let M 13 be a real, symmetric n x n matrix with non- 
negative eigenvalues. Further let Pi (A, k) be a complex n- vector, which is a function of the mode labels k and 
A, Also, let f(t) be a smooth, non-negative function on R which decays rapidly at infinity, with pointwise 
square root f l ' 2 (t) = \//(£) and Fourier transform given by 

/oo 
dtf(t)e-^. (64) 

-oo 

Then in an arbitrary quantum state |0), the following inequality holds, 

ReJ2 I d 3 kd 3 k' {/[w(k')-w(k)](4(k)av(k'))Pi(A,k)M iJ 'P i (A',k / ) ± 

±/[w(k) +w(k')] (a A (k) av (k'))P(A,k)M^P ; -(A',k')} 

>-— / */£) /d 3 k /V2[i/ + W (k)] P(A,k)AP>P,(A,k). (65) 

^ Jo x J 

The above inequality is a generalization of the scalar field positivity lemma derived by Fewster and 
colleagues |L||lE| . In order to prove this relation, first define the vector operator 

[Qt] t = E f d3k {a\ v ~ «(k)]oA(k)P,(A, k) ± g[v + w(k)]4(k)P7(A7k)} , (66) 

A J 

where 

g{u>) = -J=/V 2 M- (67) 

From the definition of the convolution 

/oo 
duj'h 1 {Lu-u')h 2 {uo'), (68) 

-OO 

it follows that {g * g) = /. 
Next, note that 

n 
a=l 

where the V?-. are the eigenvectors of M lJ , and the k q > are the corresponding eigenvalues. Now 



a=l 



ml M ij [QfU = £ « a ([Qt]\ v? a) v* a) [QtU = £ - Q l^) [OS i IV')! 2 > o • 



(70) 



Furthermore, using the commutation relations and symmetrising the integrand in (A, k) and (A',k'), we find 



J dv{[Qi]\M^ [Q±] .) = (S±) +y ^Ey rf3k lff[" + «« 2 P^k)M«P,(A,k), 



(71) 



where 



(S* ± ) = Re£ / rf 3 kd 3 k' F(k,k0(4( k )Mk0)Pi(A>k)^'P,-(A',k') 

A, A' ^ 

± G(k,k')<aA(k)a A ,(k')) P i (A,k)M i J'P i (A',k / )] 
and the functions F and G are given by 

F(k, k') = J"" dv [g[v - «(k)] .g[^-c(k')] + ffFM^I^^ + w(k'))} 



and 



The expressions for F and G may be simplified to 



G(k, k') = f°° dv [g[u + «(k))] «,[!/ - w(k')] + flf[i/ - w(k)] fl [i/ + w(k')]} 



(72) 



(73) 



(74) 



and 



F(k,k')=/ di/5[u;(k / )-i/]s[i>-u;(k)] 

•/— OO 

= (< 7 *< 7 )[u,(k')-u;(k)] 
= /[w(k')-w(k)] 

G(k,k') = /Kk)+u;(k')]. 



(75) 



(76) 



From Eq.(|70|) we know that the right hand side of Eq.(fnJ) is manifestly positive, so we conclude that (S ± ) 
obeys the following bound 



(5 ± > >- f dvJ2 I d3k |.g[^ + ^(k)]| 2 P(A,k)M y p(A,k), 



2tt 



— / dvJ2 d 3 k f^[v + Lo(k)} P(k,A)Af^P(A,k) 



(77) 



thus proving Eq. (pq). 



IV. THE QUANTUM INEQUALITY 

Consider a stationary observer whose four-velocity is given by 

« M = (lffoor 1/2 , 0,0,0). 



(78) 



In both the simple quantization scheme using the Coulomb gauge and in the Gupta-Blculcr quantization 
scheme, the energy density measured by this observer is given by the Maxwell portion of the stress-tensor, 
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p = T^ axweU u^u u 
_ 1 

~~ 2 



-l-Fiog^Fjo + ^^g^Ftr, 



Now make the identification 



Ei = F i0 and B h 



-.Qijkrjk , 



(79) 



(80) 



where Qijk is the completely antisymmetric Levi-Civita symbol. The energy density can then be written as 



P=\\9\ 



-1/2 



E t &Ej + Bi (e ij ) Bj 



where e is an ordinary 3x3 matrix with elements 

e = e(x) 



, . 9 11 9 12 9 13 

9 I ^21 ^22 ff 23 
.'/mi \ g 31 g Z2 ff 33 



(81) 



(82) 



The definitions of Ei, Bi and e have been shown to recast the curved space Maxwell field equations into 
the form of the Maxwell equations inside an anisotropic material medium in Cartesian coordinates. In this 
interpretation, e plays the role of the dielectric tensor in the constitutive relations. We will not push this 
interpretation any further and refer the reader to references |?8|-pO| for further discussion. 

Upon substitution of the mode function expansion into the stress-tensor, and making use of constitutive 
relations and the commutation relations we find 



XX' 



= \g\- 1/2 I ReJ2 / d 3 kd 3 k' \a{(k)a x ,(k') E i (A,k;x) e ij Ej(X' ,k / ;x)e 



i(u—uj')t 



-a A (k)a v (k') ^(A,k;x) & ^(A',k';x)e 



-z(cj+u/)t 



A A' 



+Re^ / d 3 kd 3 k' U(k)av(k')B,(A,k;x) (I")" 1 B^A'.k'jx) e^-"')* 
-a x QL)ax.(y)Bi(\,k;x) (e^y 1 5,(A',k';x) e -H-W)t] 



where 



and 



+ -^/<i 3 k U(A, k; x) e« EjiX, k; x) + B. t (A, k; x) (e^)- 1 Bj(X, k; x) 



^(A,k;x)=9 i [7 (A,k;x)+ia;(k)C/ i (A,k;x) 



Bi(X, k; x) = Qiji dj Ui(X, k; x) 



(83) 



(84) 



(85) 



The last line of Eq. (|83|) is the vacuum self-energy of the photons. As was the case for the scalar field, we 
will look at the difference between the energy in an arbitrary state relative to the vacuum energy using the 
normal order prescription, i.e., 



p:= p-p 



vacuum- 



(86) 



It is now our intention to show that given a temporal sampling function /(£) then the sampled energy density 
defined by 



Ap= / dt(:p(x,t):)f(t). 



(87) 
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is bounded from below. Using Eqs. ( p3[ ) and (|8q), along with the definitions of F(k, k') and G(k, k') given 
by Eqs. ([73]) and (JT^) , the sampled energy density is 



Ap = 



1 



Isool 



AA' 



ReJ2 / d 3 kd 3 k'u>(k)u;(k') |>(k, k r )(a{(k)a x ,(k')) J5<(A,k;x) £« £^(A',k';x) 



5oo 



+G(k,k')(oA(k)av(k')) ^(A,k;x) g*i ^(A',k';x) ] 
Re^2 f d 3 kd 3 k' F(k,k')(4( k ) a A'(k')) Si(A,k;x) {g lj )' 1 Bj(k',A';x) 



AA' 



+G(k,k')(a A (k)a A ,(k')) B 4 (A,k;x) (.g^)- 1 B J -(A / ,k y ;x)] . 

Clearly, both parts of the above expression are of the form (S^), so we may apply the bound (pq) with 
M ij = g 13 and P 2 (A,k) = £ 4 (A,k;x) for the first part of the expression and M ij = {g ij )~ l and P t {\,k) = 
-Bi(A,k;x) for the second part of the expression. This yields a difference inequality of 



Ap> 



2tt 



d "J2 / rf3k / l/2 [^+^( k )] 

A •* 

goo 



Isool 



B i (A,k;x) fl i«£? i (A,k;x) 



B l (A,k;x)(^")- 1 S J (A,k;x) 



(89) 



This expression is similar in form to the mode function expansion of the scalar field quantum inequality 
found by Fewster and colleagues pILO. The quantum inequality, Eq. (0), is found by adding the suitably 
renormalized vacuum energy density to the above expression. 



V. EXAMPLES 



A. Minkowski Spacetime 



This quantum inequality is easily evaluated in Minkowski spacetime with no boundaries. Using quantiza- 
tion in the Coulomb gauge, the four-vector mode functions are 



A a (A,k;x,t) = (0,A(A,k;x,t)) 



where 



A(A,k;x, t) 



A i(k-x— tot) 



v /2w(2tt) 3 



:£ t e 



(90) 



(91) 



e£ is a unit electric polarization vector and uj = Vk • k. Due to the Coulomb gauge condition, the propagation 
vector is orthogonal to the polarization vector, i.e. 



k-ii = 0. 
A third, orthogonal unit vector along the magnetic field direction is defined by 

&£ = k x g£. 
Inserting the mode functions into Eq. (p9|), and using g y = <5 y , we find 

1 r°° ^^ r . r - — - -, 2 



(92) 



(93) 



P> 



(2tt) 4 



(2tt) 3 



dvJ2 d 3 k[f/i[v + u{k)]} W (k) 



et • £k + bi ■ bi 



dv I dial 



{/Va[ v + w ]} S 



(94) 
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where we have made a change of variable in the momentum integration to spherical coordinates and have 
already carried out the angular integration and summation over polarization states. The next step is to 
make another change of variable 



v + LU, 



w, 



(95) 



to find 



P> 



(2^) 3 Jo 
1 



du 
du 



/V*( tt ) 



dv v 3 , 



V 2 /V2( W ) 



2(2tt) 3 
Using Parsaval's Identity, the quantum inequality is found to be 



P> 



dt 



d 2 



f 1/2 (t) 



(96) 



(97) 



This is the most general expression for the quantum inequality in Minkowski spacetime with an arbitrary 
sampling function. For the choice of a Lorentzian sampling function (0) with characteristic width to, it is 
straightforward to calculate 



P> 



27 



1024 ir 2 4' 



(98) 



This is a slightly stronger result, by 9/64, than the inequality proven by Ford and Roman J?J using an 
alternative method. Comparison with the quantum inequality for the scalar field in Minkowski space derived 
by Fewster and Eveson |14f| , shows that the electromagnetic field quantum inequality in Minkowski space 
always differs by a factor of two. This is a result of the electromagnetic field having two polarization degrees 
of freedom, unlike the scalar field which has only one, and the both the scalar and electromagnetic field 
modes having the same energy spectrum. Electromagnetic field quantum inequalities for various sampling 
functions are summarized in Table 1. 



B. Rindler Spacetime 

Next, we would like to find the quantum inequality in Rindler spacetime. We begin with the Minkowski 
space length element, 



ds 2 



-dt 2 + dx 2 +dy 2 + dx 2 . 



Next we apply the coordinate transformation 

t = £sinh?y, 
x = £ cosh 77, 

to arrive at the Rindler length element 

ds 2 = -i 2 drf + d£ 2 + dy 2 + dz 2 . 



(99) 



(100) 



(101) 



In this form, the metric is static, but the goo component is not a constant, so we can not quantize the 
theory in the Coulomb gauge but must use the Gupta-Bleuler formalism. Thus we are looking for mode 
solutions to the vector wave equation (B3). These, have been calculated by Candelas and Deutsch pit] for the 
two physical polarizations. The unphysical solutions have also been calculated |32|,[33[]. The modes can be 
conveniently expressed in terms of the mode solutions to the massless scalar field wave equation in Rindler 
space, 



13 



1 -32 , 1 a < : i : )2 



^ + -d^d i + d^ + d 2 z ) ( j ) (x) = r, 



(102) 



The positive frequency scalar mode solutions, normalized for the Klein-Gordon inner product of scalar fields, 
are J34| 



(ui, ky,k z ', x) 



(2tt)* 



(wahum) 1 ' 2 K iu (JlS)e« k * 



y-\-k z z—ior}) 



(103) 



>V2 



where /3 = (fc 2 + fcf) and if^a;) are the modified Bessel functions of the second kind (Macdonald func- 
tions) of imaginary order. The two physical modes that are important to our calculations are the transverse 
electric modes (TE), 



A a (l,u>,ky,k z ;x) = - (0, 0, d z , -d y )<j)(uj,ky,k z ;x), 



and the transverse magnetic modes (TM), 
A a (2, u), k y , k z ] x) = 



-q K<%, -d ni 0, 0) (j>(uj,ky,k z ;x). 



(104) 



(105) 



These two modes are properly normalized with respect to Eq. (f!2j) and are also orthogonal. If they are 
inserted into Eq. ( p9J ) for the difference inequality, and after a little algebra, we find 



-i poo POO 

A/5 > / dv I du> 

I" Jo Jo 



IT 4 , 

^ JO 



dv I duj 
lo 



/V2[l/ + W 

/V a [i/ + w 



iV* 



Qihjy Qihjg 



f44+—44 + {d i (l>){d l :<l>) 



, ., , t -^ £d ? ) 8inh(7rw) / d/3 ii Kl (00 , 



£ 2 4£ 



(106) 



where we have switched to polar coordinates to carry out the angular portion of the dk y dk z integrals. With 
the aid of Eq. 6.521.3 of pa], it is easily demonstrated that 



/■oc 

J Q 2E, Z smh(7rw) 



(107) 



Thus 



A/5>--— W dv duouj{uo 2 + \) /Va^ + o,] 
^ t, Jo Jo 



1 



16tt 3 £ 4 

'8^ < 



2 / 1/2 W 



di] 



2 f 1/2 (v) 



du + 2 



dry + 2 



*/i/a(u) 



du 



di] 



f 1/2 (v) 



di] 



(108) 



where we have again changed the variables of integration in accordance with Eq. ( |95| ) in the second line and 
used Parsaval's Identity to arrive at the third line. The quantum inequality is found by adding the Rindler 
space vacuum energy density 131 , 



1 



J vacuum 



TT 3 ^ 



dio 



Ld 3 +UJ 



11 



,2™ _ 1 240tt 2 £ 4 ' 



(109) 



to the above expression. For the Lorentzian sampling function, Eq. (0), and the definition of the proper time 
of the static observer, r = £77, we find 
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27 
Ap> 



1024tt 2 7 



27 U 



(110) 



Once again we find that the Rindler space difference inequality for the electromagnetic field is twice that 
of the scalar field result found by Fewster and Eveson |14[ for the same reason as discussed in the previous 
example. The electromagnetic field quantum inequalities for other sampling functions are also summarized 
in Table 1. 

We need to be careful about the interpretation of this quantum inequality in Rindler spacetime, as i t 



appears that both the vacuum energy density and the expression for the difference inequality, Eq. ( 108 ), 
diverge in the limit as £ — ► 0. This does not mean the quantum inequality fails on the particle horizon 
in Rindler spacetime. This divergence is really a pathology of the coordinates and spacetime trajectory 
used. Recall that the quantum inequality found above is for a static observer in the Rindler coordinates. 
This trajectory is not that of a geodesic observer but one undergoing constant acceleration. A "static" 
observer at £ = would require a constant infinite acceleration, an impossible scenario. The divergence in 
the quantum inequality expresses this. We can then ask what is the quantum inequality along the worldline 
of a geodesic observer in Rindler space? Well, a geodesic observer in Rindler spacetime is the same as a 
constant velocity geodesic in Minkowski spacetime, with the resulting quantum inequality in the geodesic 
observer's rest frame already found in the preceding Minkowski space example. It is obvious that there is 
nothing "unique" happening as the geodesic observer crosses the point is space which is associated with the 
particle horizon in Rindler coordinates. Thus, in Rindler space, the quantum inequality along a geodesic 
worldline does not fail. 



C. Static Einstein Spacetime 

Finally, we study the quantum inequality in the static closed universe where the length element is given 

by 

ds 2 = -dt 2 + a 2 [d X 2 + sin 2 X (dO 2 + sin 2 6 dip 2 )] , (111) 

and a is the radius of the universe. The modes of the electromagnetic field in this spacetime have been 
studied by various authors |3(]|57],[3(| . In terms of the vector potential, the mode solutions are the vector 
spherical harmonics on S 3 with harmonic time dependence. In a fashion similar to the previous examples, the 
four-vector mode functions can be found from a scalar function that satisfies the partial differential equation 

2 cos v \ 
V a V a -^- JL d x )^ nlm (t, x ,8,<p) = 0, (112) 

o sinx / 

which is not the scalar wave equation in the Einstein universe. The scalar mode solutions are 

ipnim(t,x,0,v) = V nl (x)Y lm (9,<p)e- i "™ t , (113) 

where u> n — n/a and Yi m (9, ip) are the scalar spherical harmonics on S 2 . The functions V n i{x) are defined as 



Vm(x) = 2 '"^ ( " } =| sinW X CW_ l( cosx), (114) 

y/l(l + ljTT(n + l)\ 

where C^(x) are the Gegenbauer polynomials as defined in B3]. The primary quantum number n ranges 

over the integers greater than one, i.e. n = 2, 3, 4, For a given n there are n 2 — 1 harmonic states with 

the same energy labeled by the quantum numbers, I = 1, . . . , n — 1 and n — —I, — I + 1, . . . , 0, . . . 1 — 1,1. 
The two physical four-vector potential modes are the electric J-pole modes, 

A a (l,n,l,m;x) = - ( 0, ——* — ,d x dg,d x d v ) il} n im(t,x,6,<p), (115) 

n \ sin x J 

and the magnetic J-pole modes, 
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A a (2,oj, k y ,k z ;x) = [0,0,- <9 v ,sinx<9e I VWi(*,X, #,¥>), 



(116) 



both of which satisfy the vector wave equation (|l0|), the Lorentz gauge and Coulomb gauge conditions, and 
are orthonormal. 

Inserting these mode into Eq. (189) yields 



Ap>— dvY, / 1/2 KwJ E 



n=1 



j-^ a 4 sin 2 x 



1(1 + 1) 



\& x + 2r? 



vti E y '™ y '" 



m=—l 



(d x V nl ) 2 + n 2 V^ 



suit 



a e sin^ + -^-a 2 ) E Y ^ Y i 



sin 



(117) 



m= — I 



However, the spherical harmonics satisfy an addition theorem, 

l 



/ j *lm*lm 



m=—l 



21 + 1 

4tt ' 



(118) 



which is independent of the 6 and ip coordinates, thus the terms in the expression for the difference inequality 
involving derivatives with respect to 9 and tp will vanish. The remaining terms can then be written as 



Ap> 



1 



27r 2 a 4 



*'£|/ 1/a ["- 



2 / 2 I 

n + 



4sin 2 x 



9>^X £ 



(2i + l)l(/ + l) 



8=1 



sin 2 x 



K 



/// ' 



(119) 



The Gegenbauer polynomials also satisfy an addition theorem, Eq. 8.934.4 of [p5[ , which for our case can be 
written as 



^ (^ + l)2 2 ^!) 2 (n-/-l)! r , , l+1 , 

2^ 77T77J \- sm xCn-i-i( C0S X)\ = n. 



1=0 



(l + l) 



(120) 



Using this in the difference inequality leads to 

Ap>--3-/ dvJ2\f 1/2 W + ^} 



2vr 3 a 3 



o 



1 



2 "X 



d^w> 2 -l) f l ^[v + LO n ] 



4 sin z x 

2 



dl sin 2 x 



1 / sin n\ 
n 



n \ smx 



(121) 



The resulting expression is spatially invariant, as expected in a homogeneous and isotropic universe. In 
addition, it has the general form of a sum over all the energies times the multiplicity for each energy times 
the Fourier transform of the square root of the sampling function, a form similar to that found by Fewster 
and Teo fig ] for the scalar field in both the three-dimensional closed universe and in the four-dimensional 
static Robertson- Walker spacetimes. In order to find the quantum inequality, we need to add to the above 
expression the renormalized vacuum energy density for the electromagnetic field which is found to be [B8fl 



11 



^vacuum — 



2407r 2 a 4 ' 
When the difference inequality is evaluated for the Lorentzian sampling function we find 



Ap> - 



27 



1024tt 2 £ 4 



SEM(to/a), 



(122) 



(123) 



where Sem( z ) is the scale function for the closed universe given by 
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Sem(z) = ^| z* V n(n 2 - 1) / X» du , (124) 

^ n=2 ^ 

and K (u) is the zero-order modified Bessel function of the second kind. It is straightforward to evaluate 
this function numerically and is plotted in Figure m. For sampling times very small compared to the radius of 
the universe, the scale function is approximately one, for which we effectively recover the flat space quantum 
inequality. This makes sense because over such sampling times the region of the universe over which the 
observer moves is indistinguishable from Minkowski space. However, for sampling times on the order of, or 
larger than the radius of the universe, the observer (and thus the quantum inequality) has time to "see" the 
large scale structure of the universe. Thus the scale function changes appreciably away from one. 

It should also be pointed out that unlike the Minkowski and Rindler spacetime examples, the quantum 
inequality for the electromagnetic field is not simply twice that of the scalar field quantum inequality. In 
both of the previous cases, the spacetimes are flat with the Riemann curvature term in the wave equation 
vanishing. Therefore, the electromagnetic wave equation (nffl) in the Lorentz gauge can be reduced to the 
scalar field wave equation. Thus, the energy spectra are identical for the scalar and electromagnetic fields 
in each spacetime with the factor of two coming from the degeneracy of electromagnetic field having two 
orthogonal polarization states. However, for the Einstein universe, and in curved spacetimes in general, 
the energy spectrum for the scalar and electromagnetic field modes are not the same, thus the scalar and 
electromagnetic quantum inequalities have different forms. 

Using the work of Fewster and Teo J15[ , the scalar field difference inequality in the Einstein closed universe 
with a Lorentzian sampling function is 

27 

A P > ~ 20 487r 2 t 4 S scalar(t /a), (125) 

where 

S SC aiar{z) = —^ z 4 V y/n(n + 2)(n + l) 2 / Ki(u) du. (126) 

2'tt ^ Jj^+T)z 

This scale function is also plotted in Figure U\ where we again see the generic behavior of the scale function 
being one for small values of to/a and decaying for large values. However, unlike the electromagnetic case 
which is monotonically decreasing function, the scalar case has a bump which peaks at to/a ~ 0.75 and then 
smoothly decays. The bump is due to the n = 1 term in the summation, a term which has no electromagnetic 
counterpart. If this term is removed from the summation, the remaining portion of the scale function does 
result in a monotonically decreasing behavior more akin to, but not exactly like the electromagnetic case. 
At present, it is not know if the bump in the scalar case has any physical meaning, as no state has yet been 
demonstrated which actually achieves this bound, although it may be a good guess that such a state would 
include n = 1 modes. There has also been an alternative conjecture that the bump may be an artifact of the 
inequalities not being optimal. In either case, further research on the scalar field quantum inequality should 
eventually clarify this issue. 
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Sampling Function 


Minkowski Spacetime 


Rindler Spacetime 


Lorentzian 
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/5^ 27 


A P - 10247r 2 r 4 


_1 + i(f) 2 ] 


n(t*+t$ 


P — 10247T 2 ** 


Lorentzian 2 


2t 3 a 


f)^ 3 


Ap> Tg^[ 


i + llf) 2 ' 




T(* 2 +*g) 2 


P — 167T 2 *4 


Gaussian 


1 e -(*Ao) 2 


f)^ 3 


A ^ 32^[ 


1 + f(f) 2 




P - 32tt 2 t i 


Cosine 4 


(4 cos4 fe) ~to<t<t 
[ elsewhere 


~P>-& 


Ap>"4 


! + ^(f) 2 " 





TABLE 1. Electromagnetic field quantum inequalities in Minkowski spacetime and difference inequalities in 
Rindler spacetime calculated for various unit area sampling functions. 
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FIG. 1. Plot of the scale functions for a Lorentzian sampling function in the four- dimensional static Einstein 
universe. The solid line is the Electromagnetic field result, while the dotted line is the scalar field result. Note, for 
small z = to/a both scale functions approach one, while for large z they decay to zero. 
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